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Abstract— The paper deals with the simulation of network 
remodeling in linear case and demonstrate how virtual 
distortion generated a chosen branch (e.g. in the branch No.4) 
can simulate the network modification due to total blocking flow 
in this branch. To this end, the condition of flow vanishing in the 
branch under remodelling should be postulated, where the 
resultant state of flow redistribution is calculated from the 
formulas superposing linear response of the original network 
configuration and the component induced by unknown virtual 
distortion. Then making use of the analytical network model [1, 
3, 4] of this installation and using presented below, the so-called 
Virtual Distortion Method (VDM), simulation of network 
remodeling can be performed. 

Index Terms — Water supply, simulation, distribution 
systems, VDM based design. 

I. Introduction 

Global demand for water is continuously increasing due to 
population growth, industrial development, and 
improvements of economic conditions, while accessible 
sources keep decreasing in number and capacity, moreover, 
the applications involving manipulation and transport of 
water and fluids in general demand high power consumption. 
The optimal use of such water supply networks seems to be 
the best solution for the present and thus it is necessary to 
carefully manage water transfer [8, 9]. 

The growing complexity of building projects and the need for 
understanding social and economic consequences of its 
conception, design, execution, and post construction usage 
call for good planning, utilization of resources, cost 
optimization and post construction performance appraisal [2]. 
Water network simulation provides a fast and efficient way of 
predicting the network behavior, calculating pipe flows, 
velocities, head-losses, pressures and heads, reservoir levels, 
reservoir inflows and outflows and operating costs [5, 6, 7]. 
The problem of the management of water resources is more 
and more important on a world scale. In particular, there is a 
requirement for novel concepts helping to solve the water 
management problem i. e. numerically efficient tools 
supporting optimal design (redesign) process for water 
networks. 

The proposed methodology is based on the so-called Virtual 
Distortion Method (VDM) approach, applicable also in the 
problem of damage identification through monitoring of 
piezo-generated elastic wave propagation [5]. This technique 
(called Piezodiagnostics) is focused on efficient numerical 
performance of inverse, non-linear, dynamic analysis. 
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The crucial point of the concept is pre-computing of structural 
responses for locally generated impulse loadings by unit 
virtual distortions (similar to local heat impulses). These 
responses stored in the so-called influence matrix allow 
composition of all possible linear combinations of the 
influence of local non-linearities (due to defect) on final 
structural response. Then, using a gradient-based 
optimization technique, the intensities of unknown, 
distributed virtual distortions (modelling local defects) can be 
tuned to minimize the distance between the computed final 
structural response and the measured one. 

II. FORMULATION OF THE SIMULATION PROBLEM 

The so called Virtual Distortion Method (VDM, used 
originally in redesign of structural systems) [4] is the basis of 
analytical formulation of the simulation problem. 

The main advantage of the proposed VDM approach is it’s 
numerical efficiency in modelling of non-linearities and 
modifications of primary design. 

For water distribution system in a steady state flow, let me 
define the following set of equations. The equilibrium of the 
system relating the internal flow distribution in the network’s 
branches Q [m /s] with the external inlet/outlet q [m /s] is 
expressed as follows: 

q = NQ (1) 

Where, 

q = external inlet/outlet flow, 

Q = internal flow in element and 
N = incidence matrix defining topology of the 
network (having only three values 0, -1, 1) 

connecting branches to nodes and showing the 
direction of flow in the network. The -1 (inlet flow) 
and 1 (outlet flow) entries are incident to the 
direction of flow. 

The relation between water head and head losses 
can be presented as follows: 

h i = N T H i (2) 

Where, h = energy loss, 

H = water head and 
N = incidence matrix. 

The constitutive relation for water networks relates the 
pressure head h with the flow Q in the elements can 
expressed as: 
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fyf (3) 

Where, Qi = flow in that element 

hi = energy loss in element i 

Ri = constant depending on pipe diameter, length, 

type. 

Substituting Eqs. (3) and (2) into (1), the following 
formula can be obtained: 

N (R N T H') 2 — q (4) 

The relation (4) is non-linear, nevertheless, let me 
temporarily assume the linearity of this relation i.e. 

NRN r H = q (5) 


m. VDM-based simulation of parameter modification 

To make the simulation of parameter modification let me 
consider an example of the water distribution network shown 
in Fig. 1. 



Fig. 1 Oriented graph modelling a 2-loop water 

network 

It has been assumed that the network is supplied 
only through the node No.l (inlet with intensity qO 
and the only outlet is through the node No.4. 

Analogously to the Virtual Distortion Method (VDM) 
applicable to the truss structures [4,5] let me postulate that 
local modification of a network parameter can be introduced 
into the system through the virtual distortion c° , incorporated 
into the formula (4): 

NR(N t H- 8°) = q (7) 

The virtual distortion 8i° is of the same character as the water 
head (see Fig. 2) and its physical meaning is an additional 
water head externally forced in branch 64 i 44 (e.g. due to a 
locally installed pump). 



Fig. 2 Distortion simulating water flow (pressure head 
modification) in branch No. 4 

The influence of virtual distortions on the resultant flow 
redistribution can be calculated using the so-called influence 
matrix Dy collecting i responses (row-wise) in terms of water 

heads ~ l induced in the network by imposing the unit 
virtual distortion 

£ °j=l generated consecutively in each network branch j. 

Thus each influence vector H® 6 -1 can be calculated on the 
basis of the following equation obtained from Eq. (7): 

NRN T tf ' = q* + NRI (8) 

The vector q* disregards the external inlet and outlet (the flow 
is now provided by the imposition of virtual distortion), and it 
accounts for the water flow distribution in the closed network 
(cf. Eq. (6)). There is a set of j (j the number of branches) 
equations (8) to be solved in order to create the full influence 
matrix D. Each time the right hand-side changes as the unit 
virtual distortion is applied to another branch. In practice this 

can be realised by applying a pair of inlets-outlets L ik R k j £°j 

corresponding to each branch (cf. Eq. (7)) - it is the so-called 
compensative charge. 

So, the parameter modification in the system is accounted for 
by superposing the so-called linear response of the original 
network and the so-called residual response due to imposition 
of the virtual distortion. Therefore, the resultant water head 
distribution can be expressed as: 

H,=H, L +Hf =Hj-+ID^ (9) 

j 

And the resultant water flow as: 

Q J =Q^ + Q?=Qj+R l m(D i| -a ik do) 

j 

The analogous set of relations governs the VDM based 
approach to modifications of truss structure system [3]. 
Coming back to the example shown in Fig. 2, let us generate 
the unit virtual distortion in branch No. 4, connecting the 
nodes Nos. 2 & 3. The corresponding set of equations (8), 
accounting for boundary conditions (i.e. outlet in node No.4), 
takes the following form: 
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where 84 = 1. Assuming the following data: K^O.2 m 3 /s, 

K 2 =K 3 =K 4 =K 5 = 0.4 m 3 /s, 1 1 =1 2 =1 3 =1 5 =10.000 m, 1,=14.142 m, 
qi=0.050 nr/s, H 0 =0.000 m, we get the following set of 
equations for the water head distribution: 
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The resulting distribution of water heads H @s -1 = [0.151, 
-0.251, 0.251, 0.000] T constitutes the 4th column of the 
influence matrix D. Continuing this procedure for virtual 
distortions generated in other branches, the full influence 
matrix can be determined as: 
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( 12 ) 


Taking into account relation (3) and applying it consecutively 
to each influence vector H @e° 

1 , another influence matrix D £ 
can be created, collecting the response to unit virtual 

distortions in terms of the pressure head 8 @ £ _1 . 
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(13) 


Si = sf' + Z Dij s°j 

, 1 (14) 

Q=Q, +R,T(iy;-8, j )s] 

j 

Therefore, the virtual distortion to be generated in branch 
No.4 to simulate complete blocking of local flow can be 
calculated from the following condition: 

Q 4 = Q 4 + Ri (Du - 1) S 4 = 0 , or making use of (4) 
£4 ' Dij £4 £4 


What leads to: 



= 1,34 m 


(15) 


Finally the pressure head as well as the flow in modified 
network is (after substitution value (15) to relations ( 6 )) as the 
following: 


si = s[ + DU si = 3,04-396* 134 = 3,57m 
S 2 = sh + DUsl = 2,365-099*1,34 = 2,23m 
ft = ft + Dm ft = 1,225 - 0,247 * 1,34 = 0,89m 
ft = ft + DU ft = 1,9 + 0,248* 1,34 = 2,23m 

and the flows: 

Qi = +/?, Z)[ 4 ft = 0,01216+ 0,004 * 0,396 * 1,34 = 0,0143 m 3 /s 
Q 2 = + ft DUsi = 0,03784+ 0,016(-0,099)1,34 = 0,0357 

m 3 /s 

Q :i = Q \+/? 3 Df 4 ft = 0.0196+ 0.016(^0.247)134 = 0,0143 m 3 /s 

Q. = Q \+ft UiJl = 0.0304+ 0.0 1 6x 0.248x 1 .34 =0,0357 m 3 /s 

For comparison, let me solve the set of equations (14) 
taking into consideration excluding the element No. 4 (i.e. 
assuming R 4 = 0 and disregarding column 4 in the matrix L) 
one can get the following set of equations: 


IV. SIMULATION OF LINEAR REMODELING 
(ELEMINATION OF BRANCH 4) 

First, let me demonstrate how virtual distortion 
generated a chosen branch (e.g. in the branch No.4) can 
simulate the network modification due to total blocking flow 
in this branch. To this end, the condition of flow vanishing in 
the branch under remodelling (Q 4 = 0) should be postulated, 
where resultant state of flow redistribution is calculated from 
the formulas superposing linear response of the original 
network configuration and the component induced by 
unknown virtual distortion: 
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The resulting distribution of water head is: H = [4.514 
0.943 2.282 0.05] T , which leads to the following state of 
pressure head as well as the flow in modified network (after 
substitution H’ to (3) and (14)): 
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81 = Hi - H 2 = 4,514- 0,943 = 3,57m 
g2 = ni - H 3 = 4,514- 2,282 = 2,23m 
83 = H 2 - H 4 = 0,943- 0,05 = 0,89m 
g5 = ni -H 4 = 2,282-0,05 = 2,23m 

and the flows: 

Qy = Qy+Ri D £ U a = 0,01216+ 0,004x0,396x1,34 =0,0143 
m 3 /s 

Q 2 = Q 2 +R 2 Du£ ( l=Q, 03784+0,016(-0,099)xl,34 = 0,0357 m 3 /s 
Qy= Q\+Ri Dm si= 0,0196+0,016x (-0,247) x 1,34=0,0143 m 3 /s 

2 5 = Q\+Ri DUsl = 0,0304+ 0,016x 0,248x 1,34 = 0,0357 m 3 /s 

States (3) as well as (4) are the same, what demonstrates that 
virtual distortion (14) models properly the assumed 
modification. 

Multiplying the network response H R =[-0.151 0.251 -0.251 
0.000] for the unit virtual distortion ^ = 1 7 by the 

determined above value, s° 4 the searched correction to 

pressure distribution in the original network can be calculated 
(Fig.3c) and the resultant pressure distribution for the 
modified network: H = H L + 1.34 H R can be also determined 
(Fig.3b). Similarly, other, various types of the network 
modifications can be simulated through virtual distortions 
using determined once the initial system matrix, the linear 
response H L and the influence matrix D £ . 



Linear constitutive relations have been assumed in the above 
considerations. In the case of nonlinear problem formulation a 
superposition of two virtual distortion fields has to be taken 
into consideration. The first one, 8° modeling system redesign 
and the second, p° modeling physical nonlinearity of the 
system. 
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Fig.3 Pressure distributions for the original (a), locally 
distorted (b) and modified (c) networks 

For large water networks with small, local modifications the 
above VDM based approach is much cheaper numerically 
than the classical way through recomposing and solving the 
modified system. 


VI CONCLUSION 

The main advantages of the VDM based approach to the 
water network analysis are the reduction of numerical costs 
and avoidance of iteration due to incremental approach in the 
analysis of water network. 

The numerical cost of linear analysis consists of: 

Solving the linear problem (7). 

Composing the influence matrix D s (13). 
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